Light Stop from b-r Yukawa Unification 
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We show that b-r Yukawa unification can be successfully implemented in the constrained minimal supersym- 
metric model and it yields the stop co-annihilation scenario. The lightest supersymmetric particle is a bino-like 
dark matter neutralino, which is accompanied by a 10-20% heavier stop of mass ~ 100 — 330 GeV. We high- 
light some benchmark points which show a gluino with mass ~ 0.6 — 1.7 TeV, while the first two family 
squarks and all sleptons have masses in the multi-TeV range. 
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The apparent unification at Mq w 2 x 10^'' GeV of the 
three Standard Model (SM) gauge couplings, assuming TeV 
scale supersymmetry (SUSY), strongly suggests the existence 
of an underlying grand unified theory with a single coupling 
constant. The minimal supersymmetric SU(5) model, in addi- 
tion to unifying the gauge couplings, also predicts unification 
at Mq of the third family bottom (6) quark and tau lepton 
(r) Yukawa couplings [1]. This h-r Yukawa unification (YU) 
is to be contrasted with the minimal supersymmetric 50(10) 
and S'C/(4)c x SU{2)l x SU{2)r models which predict t-b-r 
YU [2], where t denotes the top quark. The low energy im- 
plications of b-T [3] and t-b-r [4] YU have been discussed in 
the recent literature. For instance, t-b-r YU is not realized in 
the mSUGR A/constrained minimal supersymmetric standard 
model (CMSSM) [5] because of the difficulty of implement- 
ing radiative electroweak symmetry breaking (REWSB) [6]. 

In this paper we explore the low energy consequences of 
implementing b-r YU in the CMSSM framework. REWSB 
in this case is not an issue anymore. We refer to this combi- 
nation of b-r YU and CMSSM as YCMSSM, the 'Yukawa' 
constrained version of CMSSM. Among other things, we re- 
quire that YCMSSM delivers a viable cold dark matter (DM) 
candidate (lightest stable neutralino) whose relic energy den- 
sity is compatible with the WMAP measured value [7]. One of 
our main observations is that the allowed fundamental param- 
eter space of CMSSM is strikingly reduced in the YCMSSM 
setup. We find that A/1/2 ^ ™o, where M1/2 and toq denote 
universal gaugino and scalar soft SUSY breaking masses re- 
spectively. Furthermore, b-r YU at the level of 10% or better 
yields the constraint 5 TeV < toq ^ 20 TeV. The supersym- 
metric threshold corrections including finite loop corrections 
to the b quark mass play an essential role here [8]. 

The lightest supersymmetric particle (LSP) neutralino is es- 
sentially a bino, the spin 1/2 supersymmetric partner of the 
U{1)y gauge boson, which is closely followed in mass by 
a slightly heavier (next to lightest sparticle for short NLSP) 
stop, a scalar partner of the top quark. The desired LSP relic 
abundance is achieved via neutralino stop co-annihilation, 
which in our case requires that the NLSP stop is about 10- 
20% heavier than the neutralino. The parameter tan/3, the 
ratio of the up and down Higgs VEVs, turns out to lie in a 
narrow range 35 < tan /3 < 40. The universal trilinear scalar 



coupling (Ao) is found to satisfy |^o/"^o| ~ 2.3. 

We highlight some LHC testable benchmark points with 
comparable LSP neutralino and light stop masses of around 
100 — 330 GeV, while the corresponding chargino and sec- 
ond neutralino masses are 200 — 600 GeV and gluino mass 

0.6 - 1.7 TeV. Together with the lightest SM-like Higgs 
with mass 114 - 124 GeV, these are the only 'light' (LHC ac- 
cessible) particles predicted in this NLSP stop scenario with 
b-r Yukawa unification and neutralino DM. The squarks of 
the first two families, the heavy stop, the two sbottom par- 
ticles, and the charged sleptons, all have large (multi-TeV) 
masses. 

The fundamental parameters of CMSSM are: 

mo, Ml/2, tan 13, Ao,sgn{fi) (1) 

where sgn(/i) is the sign of supersymmetric bilinear Higgs 
parameter All mass parameters are specified at Mq. 

We use the ISAJET 7.80 package [9] to perform random 
scans over the CMSSM parameter space. ISAJET employs 
two-loop renormalization group equations (RGEs) and defines 
Mq to be the scale at which gi = g2- This is more than 
adequate as a few percent deviation from the exact unification 
condition 53 = gi = 52 can be assigned to unknown GUT 
scale threshold corrections [10]. The random scans cover the 
following parameter range: 

< mo < 25 TeV, < M1/2 < 2 TeV, 

1.1 < tan/3 < 60, -3 < Aq/too < 3, (2) 

with II > and mt — 173.3 GeV [11]. The results are not 
too sensitive to one or two sigma variation in the value of nit 
[12]. 

In scanning the parameter space, we employ the 
Metropolis-Hastings algorithm as described in [13]. All of 
the collected data points satisfy the requirement of REWSB, 
with the neutralino in each case being the LSP. We direct the 
Metropolis-Hastings algorithm to search for solutions with 
10% or better b-r Yukawa unification (YU). After collecting 
the data, we impose the experimental mass bounds on all par- 
ticles [14], and use the IsaTools package [15] to implement 
the following phenomenological constraints: 
As far as muon anomalous magnetic moment is concerned, we 
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FIG. 1 : Evolution of bottom (green) and r (red) Yukawa couplings 
without (a) and with (b) finite SUSY threshold corrections. 

m/i (lightest Higgs mass) > 114.4 GeV [16] 
BR{B, n+fj.-) < 5.8 X 10-* [17] 
2.85 X 10-4 < BR{b sj) < 4.24 x lO""' (2ct) [18] 
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only require that the model does no worse than the standard 
model (SM). In Fig. 1(a) we show the evolution of b and r 
Yukawa couplings without the SUSY threshold corrections for 
a representative b-r YU solution. It is evident that without 
suitable the SUSY threshold corrections, b-r YU occurs at 
around 10* GeV. In Fig. 1(b) we show the need for SUSY 
threshold corrections in order to achieve b-r YU. 

The SUSY correction Snir to the t lepton mass is given 
by SrUr — v cos pSyr- The finite and logarithmic correc- 
tions [20] to the T lepton mass are typically small, so that 
the value of ijr at Mq is more or less fixed. From Fig. 1(a), to 
implement Ut- ~ yi, at Mq therefore requires suitable thresh- 
old correction {5yb) to the bottom quark [8]. The dominant 
contribution to Syi, comes from the gluino and chargino loops 
[8, 20]. With our sign conventions ( the sign of Syi, is fixed 
as yb evolves from Mq to Mz), a useful approximate formula 
for the finite one loop correction to 5yb is given by 



5vl 



47r2 



53^ 

3 m\ 



+ 



8 ml 



tan /?. 



(3) 



Here is the strong gauge coupling, mg is the gluino mass. 
At is the stop trilinear coupling, and toi « (mj^^ + m^^^)/2, 

7712 ~ ("T'tj + bi, 62 denote the two bottom squarks, 

t2 is the heavier stop, and we assume that 771^ <C 7n^^ and 

In order to achieve b-r YU, j/h must receive a negative con- 



tribution (—0.2 < Syt/yt 



< 



-0.07) from threshold correc- 



tions. This is a relatively narrow interval compared to the full 
range —0.2 < Syb/yi, < 0.25 in the data that we have col- 
lected. The logarithmic corrections to yb are in fact, positive, 
which leaves the finite corrections to provide for the correct 
6yii. Since fi > 0, the gluino contribution is positive, and so 
the contribution from the chargino loop not only has to cancel 
out the contributions from the gluino loop and the logarithmic 
correction, it also must provide the correct overall (negative) 



contribution to Sy^. This is achieved only for suitable large 
777,0 values and large negative At, for which the gluino con- 
tribution scales as Mi/2/mQ while the chargino contribution 
scales as At/mQ. Note that large values of 7719 imply heavy 
slepton masses, and so b-r YU does not provide any signif- 
icant SUSY contributions to the muon anomalous magnetic 
moment. 

In order to quantify b-r YU, we define the quantity i?bT as 



max{yb,yr) 
mm{yb,yr) 



(4) 



In Fig. 2 we present our results in the Rbr - mo, Rbr - tan /?, 
Ml/2 - n^o and Rbr - /i planes. In panels (a, b, d) the 
gray points are consistent with REWSB and xl LSP. The or- 
ange points satisfy the particle mass bounds, constraints from 
BR{Bs fJ-'^fJ-^), BR{Bu TVr) and BR{b 57). 
The blue points form a subset of orange points that satisfies 
WMAP bounds on x\ DM abundance. The dashed (red) line 
represents 10% or better b-r YU. In the A/1/2 - 'n^o plane the 
orange color has the same meaning as described earlier, red 
color represents solutions with 10% or better b-T YU, and the 
blue points satisfy in addition the WMAP bounds. In the Rbr 
- mo plane of Fig. 2, we show that in order to have 10% or 
better b-r YU consistent with the collider bounds, we require 
7770 ^ 5 TeV. But if b-r YU is to be compatible with neu- 
tralino DM relic density, represented by the blue points in the 
figure, we can see that mo > 8 TeV. To better appreciate this 




FIG. 2: Plots in Rbr - mo, Rbr - tan/?, Afi/2 - mo and Rbr - fJ, 
planes. In panels (a, b, d) the gray points are consistent with REWSB 
and xi LSP. The orange points satisfy in addition the particle mass 
bounds, constraints from BR(Bs — > fi^ ^i^), BR{Bu — >■ rUr) and 
BR(b — >■ 57). The blue points form a subset of orange points that 
satisfies WMAP bounds on Xi DM abundance. The dashed (red) 
line represents 10% or better b-r YU. In il/1/2 - mo plane the orange 
color has the same meaning as described above, red color represents 
the solution with 10% or better b-r YU, and blue points satisfy the 
WMAP bounds. 
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result we consider the M1/2 — rriQ plane. As we will see later, 
the neutralino-stop co-annihilation channel is the only solu- 
tion we have found for neutralino DM compatible with h-r 
YU. The lighter stop mass can be as light as 100 GcV [14], 
which means that the neutralino mass has to be of comparable 
order, in order to implement the stop co-annihilation solution. 
An ISAJET two loops analysis yields M1/2 > 150 GeV. We 
can see from Fig. 2(c) that A/1/2 ^ 150 GcV is compatible 
with 10% or better YU (points in red) for mo > 8 TeV. 

In the Rhr - tan/3 plane (Fig. 2(b)) we have two regions 
with 10% or better h-r YU. In one region tan (3 k, 1.1, but this 
is ruled out by the light CP-even higgs mass bound. The sec- 
ond region with successful h-T YU occurs for 34 < tan /3 < 
45. Beyond this range we do not have REWSB if we require 
10% or better h-r YU. Imposing , in addition, the neutralino 
DM requirement (represented by blue points), the allowed re- 
gion shrinks to 35 < tan /3 < 40. 

The plot in R}jt - fi plane shows that for 10% or better YU, 
/i > 3 TeV (orange points). The need for large values of /i can 
be understood by analyzing the threshold effects associated 
with the b quark. As previously mentioned, the second term 
in Eq. (3) has to be larger than the first one for successful 
b-T YU. With tan /3 squeezed in a relatively narrow interval, 
and the magnitude of At bounded to avoid color and charge 
breaking, needs up with a large fj, term to obtain the required 
6yii. On the other hand, a bino-Higgsino mixed DM scenario 
requires m^o « From Fig. 2(c), the neutralino cannot be 
much heavier than 600 GeV or so for successful b-r YU. But 
with /i > 4 TeV that a bino-Higgsino mixed DM scenario is 
not realized here. 

Next let us consider the A-funnel scenario of CMSSM 
where one needs ruA ~ 2m^a. For large tan/3, = 
26/sin(2;9) sa 6tan/3, where b is the soft SUSY breaking 
Higgs bilinear term. We have b tan /? « 4 TeV for good b- 
T unification, which implies that niA ^ 2to-^o and thus, the 
A-funnel scenario is not consistent with 10% or better b-r uni- 
fication. 

The neutralino-stop co-annihilation channel is compatible 
with 10% or better b-r unification, as shown by the blue points 
in Fig. 2. Let see how a light stop mass (m^^) of O(IOO) GeV, 
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FIG. 3: Plots in mi^/rrijj^ — Ao and mj^/m^^ — At planes. 
The orange points satisfy the particle mass bounds, constraints from 
BR{Bs /i^M"). BR{B^ -5> TVr) and BR{b 37). Red 
color represents the solution with 10% or better b-T YU. Blue color 
represents the neutralino-stop co-annihilation solutions. 
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TABLE I: Point 1 has perfect b-T unification with Ri,t — 1, point 2 
represents a solution with relatively light stop mass (~ 114 GeV), 
and point 3 represents a solution with a heavier stop mass (~ 
330 GeV). The remaining squarks and sleptons all have masses in 
the multi-TeV range. 

is realized for mo > 8 TeV and M1/2 < 600 GeV. In this re- 
gion of the parameter space the diagonal entries of the stop 
mass matrix to pick up dominant contribution from yt and At 
couplings, thus making the m^^ entry lighter than its value 
at Mq. In Fig. 3(a) we see that the ratio mf^/mo can be as 
small as 0.05 for large values of Aq, which means that at the 
SUSY scale m^^ can be as light as 1 TeV or so, despite the 
large ~ (8 — 20) TeV tuq values. On the other hand, from 
Fig. 3(b) we see that at SUSY scale, the value of \ At \ can be 
0(7 - 20) TeV. The off -diagonal entries ~mt{At + fi cot /?) 
for the stop quark mass matrix can therefore be of compa- 
rable magnitude to the diagonal entries. Because of this, as 
seen in Fig. 3(b), one of the eigenvalues (m^^) of the stop 
quark mass matrix can be 10-20 times smaller than m^^ at 
the SUSY scale. Thus, with m^^ O(IOO) GeV, we can have 
neutralino-stop co-annihilation compatible with b-r YU. A 
similar discussion for the stau mass matrix shows that the stau 
co-annihilation scenario is not realized in YCMSSM. 

In Table I we present three characteristic benchmark points 
which satisfy all, especially the dark matter, constraints. Point 
1 displays essentially perfect b-T YU solution with Ri,t- = 1, 
point 2 represents a solution with a relatively light stop mass 
(~ 114 GeV), and finally point 3 represents a solution with a 
heavier stop mass (^ 330 GeV). We note that the remaining 
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squarks and sleptons all have masses in the multi-TeV range. 

Since the LSP is essentially a pure bino, both its spin- 
independent and spin-dependent cross sections on nucleons 
are rather small [21], ^ 10"^'' — 10~'**cm^. Consequently it 
wont be easy to detect the LSP in direct and indirect experi- 
ments. 

In summary, we have investigated b-r Yukawa unification 
in the mSUGRA/CMSSM framework and find that it is con- 
sistent with the NLSP stop scenario and yields the desired 
LSP neutralino relic abundance. This YCMSSM predicts that 
there are just two 'light' (LHC accessible) colored sparticles, 
namely the NLSP stop with mass 100 - 330 GeV, and 
the gluino which is ^ 600 — 1700 GeV. The chargino and 
a second neutralino are about a factor 2-3 lighter than the 
gluino. The remaining squarks as well as all sleptons have 
masses in the multi-TeV range. Regarding the fundamental 
CMSSM parameters, we find that 5 TcV < mo < 20 TcV, 
mo/Mi/2 « 30 - 50, tan/3 « 35 - 40, - 3 - 15 TcV 
and lAo/mol 2.2 - 2.4. 
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